arXiv:gr-qc/0503095v2 20 Jul 2005 


Conformally coupled scalar solitons and black holes with negative cosmological 

constant 

Eugen Radi0 

Department of Mathematical Physics, National University of Ireland, Maynooth, Ireland. 

Elizabeth WinstanlejQ 

Department of Applied Mathematics, The University of Sheffield, 

Hicks Building, Hounsfield Road, Sheffield, S3 7 RH, United Kingdom. 

(Dated: February 4, 2008) 

We present arguments for the existence of both globally regular and black hole solutions of the 
Einstein equations with a conformally coupled scalar field, in the presence of a negative cosmological 
constant, for space-time dimensions greater than or equal to four. These configurations approach 
asymptotically anti-de Sitter spacetime and are indexed by the central value of the scalar field. 

We also study the stability of these solutions, and show that, at least for all the solutions studied 
numerically, they are linearly stable. 
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I. INTRODUCTION 

Within the various theories violating the no-hair conjecture, the case of a conformally coupled scalar held is of 
particular interest. In asymptotically hat space with no scalar self-interaction potential, this theory admits an exact, 
closed form black hole solution, which has a scalar held diverging on the event horizon |lj. However, it is known that 
this solution is unstable |§], and cannot be considered as a valid example of a black hole with scalar hair. There are 
also a number of theoretical results as well as numerical evidence against the existence of black holes with scalar held 
hair (with various couplings to the Ricci scalar curvature) in asymptotically hat spacetime (see, for example, |J for 
a recent discussion). 

In an unexpected development, hairy black hole solutions have been found in both theories with minimally as well 
as nonminimally coupled scalar helds by considering asymptotically anti-de Sitter (AAdS) boundary conditions 111 
Moreover, some of these solutions are found to be stable. Exact four-dimensional black hole solutions of gravity with 
a minimally coupled self-interacting scalar held have been presented by Martinez, Troncoso and Zanelli (MTZ) 1 
and Zloshchastiev Q. 

However, in many other theories admitting hairy black hole solutions, these configurations survive in the limit of 
zero event horizon radius, yielding particle-like, globally regular conhgurations. Motivated by this observation, we 
consider in this paper the case of a conformally coupled scalar held in an n—dimensional AAdS spacetime (with n > 4) 
and look for both globally regular and black hole solutions. In the black hole case, we extend the results of Ref. Uby 
considering higher dimensional conhgurations. Since a negative cosmological constant allows for the existence of black 
holes whose horizon has nontrivial topology, we consider, apart from spherically symmetric solutions, topological black 
holes also. We hnd that the spherically symmetric solutions admit a nontrivial regular limit, representing gravitating 
scalar solitons. These conhgurations are indexed by the central value of the scalar held and are found to be stable 
against linear huctuations. 

The outline of this paper is as follows: in Section El we introduce our model, the numerical results being presented 
in Section mu The stability of our solutions is addressed in Section llVl Our conclusions are presented in Section El 
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II. THE MODEL 


A. The ansatz and field equations 


We consider the following action, which describes a self-interacting scalar field cf> with non-minimal coupling to 
gravity in n— dimensions (throughout this paper we will use units in which c = 87 tG = 1) 


S = 


J d n x \f~-g 


X -[R-2A)~\ (V0) 2 - \jRcf 2 - V{<t>) 


(1) 


where R is the Ricci scalar curvature, A = — (n — 2)(n — l)/2£ 2 is the cosmological constant and £ is the coupling 
constant. For a minimally coupled scalar field £ = 0 and for conformal coupling (which is the focus of this paper) 
? = £c = (n — 2)/4(n — 1). 

The field equations are obtained by varying the action m with respect to field variables g^ and <f> 


(1 - & 2 ) G M „ + g^A = (1 - 2£) + (0 - 0 g„ v 

-2£<j>VpV u tf> + 2^ A(J/ ^V 2 ^- g^uV{<t>), 

,, dV 

V 2 


„ dV 

7*</>-ZR4-— = 0. 

The Ricci scalar expression which follows from these equations is 

2(n - 1)(£ - Cc)(V^) 2 + 2£(n - - n(V(<l>) + A) 


R = -- 


n/2 - 1 + 2 (n - !)£(£ - £ c )0 2 


( 2 ) 

( 3 ) 


( 4 ) 


Since for a negative cosmological constant topological black holes may appear (with a nonspherical topology of the 
event horizon), we consider a general metric ansatz 


ds 2 = 


dr 2 

W: 


+ r 2 da 2 -H(r)e 2S ^dt 2 , 


( 5 ) 


where dc t 2 2 k = di/) 2 + f 2 {il))dVl\ denotes the line element of an (n — 2)-dimensional space S^ with constant curvature. 
The discrete parameter k takes the values 1, 0 and —1 and implies the following form of the function 


sin if), for k = 1 

fk(4>) = { Vh for k = 0 

sinh'f/;, for k = —1. 


( 6 ) 


When k = 1, the hypersurface £i represents a sphere, for k = — 1, it is a negative constant curvature space and it could 
be a closed hypersurface with arbitrarily high genus under appropriate identifications. For k = 0, the hypersurface 
So is a (n — 2)-dimensional Euclidean space (see e.g. the discussion in J|). 

A convenient parametrization of the metric function H is 


H(r) = k-™4 

v ' r n—3 


r 

p' 


( 7 ) 


Although a rigorous computation of the solutions’ mass and action is a nontrivial task for £ ^ 0, we assume that the 
asymptotic value of m(r) corresponds to the mass of our configurations, up to some n— dependent factor. 

Within this ansatz, we find the resulting equations (where a prime denotes the derivative with respect to r) 


n — 2 
2 r 




H' - - —-( k-H) 


- ( 2£ - l ) H(f)' 2 + £#'(#' + 2H5') - 2 i(j)A7 2 ^ + + A = 0, 


n — 2 


(1 - ^ 2 ) 6 ' - (1 - 2^' 2 - 2 ^( 5 '</>' - </>") = 0 , 


H<t>" + $ ( H6' + H' + H ——- £,R(j) - = 0. 

r J d<p 


( 8 ) 
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B. Boundary conditions 


The asymptotic solutions to these equations can be systematically constructed in both regions, near the origin (or 
event horizon) and for r 1. 

The corresponding expansion as r —> 0 is (globally regular solutions may exist for k = 1 only): 


H(r) = 1 + 
<5(r) = <5 0 + 
<t>(r) = <f> o + 


2(-2em 0 + A + V 0 - 2^ 0 V 0 ') r 2 + Q ^ 


(n- l)(n- - 1) 

HMMoRo + K) 2 + Q(r 4)^ 


(n~ 1)(«- 2)(^o - !) 
^o-Ro + K r 2 + 0 


where 


i?n — 


2(n — 1) 

n(V 0 + A) - 2g(n - 1)^^ 
n/2 - 1 + 2(n - !)£(£ - £ c )</>o 


(9) 


( 10 ) 


is the Ricci scalar evaluated at the origin and Vq = V(4>o), Vq = dV/dtj) 

For black hole configurations with a regular, nonextremal event horizon at r = rh, the expression near the event 
horizon is 


where 


H\r h ) 
<t>'(r h ) = 
4*2 = 

S'(rh) = 


H(r) = H'(r h )(r-r h ) + 0(r-r h ) 2 , 

S(r) = Sh + S'(rh)(r - rh) + O (r ~ rn) 2 , 

</>(r) = (j> h + <t>'(r h )(r - r h ) + 4> 2 (r - r h ) 2 + O (r - r h f , 


, k 2r h (A + V h -£(f>h(£<t , h Rh + V/ i )) 
(n “ + - (»-2)««-D- 


( 11 ) 


ZRhfa + Vh 
H'(r h ) ’ 
4>'(r h ) 


2r h 

+£,4>h4>' (rh) 


2R;,;o,. 

k(n — 3) 


1 - 


[ r /i(£R/j + 1/T) + 2fc(n — 3)] 
rn 


m -!) 


(2g - l)^'^) 2 ^^ 

2(n - 2)H'(r fc ) 


2r h H'(r h ) 2(n - 2)ff'(r h ) 
r/t[(2^ - l)<^'(rfr) 2 + 4g^(/> 2 ] 

(n - 2)(£0? - 1) + 2 £<j) h <j>'(r h )r h ’ 


(3£ifc + o 


and 


Rft = 


n(lfr + A) - 2g(n - l)<faV/ 
n/2 - 1 + 2(n - 1)£(£ - £ c )(^ 

(k) 


( 12 ) 


(13) 


is the Ricci scalar evaluated at the event horizon; noting that V/ = (4 ) )\ < p =( p h ■ 

To analyze the r 1 region, we assume that the geometry approaches asymptotically the AdS spacetime and that 
the function m(r) does not diverge in the same limit. These assumptions imply that linv^oo <j>(r) = 0 and 


Cl 


C2 


C3 


j.n/2—l j.n/2 j.n/2+1 


H (r) - k 2M 2Af2 

V) r n ~ 3 (n — 2)(n — 1) 


+ ..., S(r) = ~A + .. 


(14) 


and also impose some constraints on the scalar potential (for example, in n = 4, the potential should satisfy as r —> oo 
the condition V = V' = V" = V"'" = 0). In the above relations ci, c 2 and M are real constants which fix the values 
of the other coefficients in the asymptotic expansion. In the simplest case of a vanishing self-interaction potential, we 
find 

cik(n — l)(n — 2 ) 2 (n — 4) 


C 3 = -- 


So = - 


16A 

c\(n — l)(n — 2) 3 (n — 4 )k + 8c|nA 
16n(n — l)(n — 2)A 


(15) 
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C. Conformal transformation 


The conformal transformation m 


g^ 2) gHU, with 12 = 1 - ^fa 2 (16) 

maps the original system 0 onto a much simpler one involving just a minimally coupled scalar field, but with a more 
complicated potential (this transformation is valid only for those solutions with a nonvanishing 12). For a conformally 
coupled scalar field, the new action principle takes the form 


S= d n x y/—g 


R — 2A 1 .2 ... . 

— 2 - 2 ^ 


(17) 


where a bar denotes quantities calculated using the transformed metric g and we define a new scalar field <1> as (for 

£ = £c): 


$ = —^arctanhy 7 !^, 


(18) 


with a nonvanishing potential 


um 


[A + V(fa)\ cosh(- v /^$) 


-A. 


(19) 


The main advantage of the rescaled frame is that the Held equations are much simpler. However, the potential m 
is unphysical (for example, with V = 0, it is negative everywhere). 

In the transformed frame, the metric g^u takes the form ©, but with the quantities H and 6 replaced by H and S 
respectively. In addition, there is a new radial co-ordinate f = f2 1 /( n-2 V, which is a good co-ordinate as long as 

A=(n- 2)12 - 2r£ c #' > 0. (20) 


This is an additional constraint on the scalar field <p required for the conformal transformation to be valid. For all 
our numerical solutions, the conditions 12, A > 0 are satisfied. 


III. NUMERICAL RESULTS 

As analytic solutions to the coupled nonlinear equations @ appear to be intractable for every dimension, except 
for the n = 4, k = — 1 black hole solution found in [f|, the resulting system has to be solved numerically. 

In this section we discuss mainly the case of a conformally coupled scalar field without a self-interaction potential. 
Since for V(fa) = 0 the field equations are invariant under the transformation </> —> —(f>, only positive values of fa 
are considered. Here fa denotes the initial value of the scalar field at the r = ro (with tq = (0, rfa) for regular 
and black hole solutions, respectively). Also, by rescaling the radial coordinate (together with m(r)), we can set 
A = — (n — l)(n — 2)/2 (i.e. i = 1) without any loss of generality. 

We follow the usual approach and, by using a standard ordinary differential equation solver, we evaluate the initial 
conditions at r = ro + 10” 6 for global tolerance 10~ 12 , and integrate towards r —> oo. In this way we find that 
nontrivial solutions may exist in any dimension n > 4 (both black hole and regular solutions exist also in three 
spacetime dimensions; however, their properties are somewhat special and we do not discuss them here). Black hole 
solutions seem to exist for any values of the parameters (fc, 4>h, rfa) satisfying H'fafa > 0. 

Typical profiles are presented in Figure 1 for regular configurations and in Figures 2, 3 for black hole solutions. 
The dependence of the mass parameter M and e 2<5 ( r °) on fa is plotted in Figure 4-6 (note the occurrence of negative 
values of M for k = — 1 black holes, a common situation in topological black hole physics). 

The properties of the configurations can be summarized as follows: 

1. For n > 4, AAdS solutions exist for any values of fa in the interval 0 < fa < 1/v^; 

2. The scalar field interpolates monotonically between fa and zero and has no nodes; 

3. The value of the metric function e 2<5 at the origin (event horizon respectively) decreases for an increasing fa and 
approaches zero in the limit fa —> 1/VCc- 
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As seen in Figure 4, for regular solutions M approaches a finite value in the same limit. In the black hole case, the 
critical value of M increases very rapidly with r^, which makes its accurate determination a difficult task for large r^. 
In the k = — 1 case, the condition H'(rh) > 0 implies the existence of a minimal value of the event horizon radius, 


( n ~~ l)( n — 3)(1 — 

[1 + n(n — 2)(1 — 0|/8)]i/2' 


( 21 ) 


We have also found that non-trivial configurations may exist in the presence of a nonzero scalar potential. In this 
case the scalar field equation © in a fixed AdS background has two exact solutions 


+ -p 


v{4>) = ccj) s , 


( 22 ) 


withp = (2 —n)/4, c = — (n— 2) 3 /8n£ 2 , s = 2n/(n—2) in one case, and p = — n/4, c =—n 3 /8£ 2 (n+2), s = 2{n+2)/n 
in the other. 

However, we have restricted our analysis to the particular form V((f>) = l/2fi 2 (f> n , which for n = 4 corresponds to 
the case considered in @j. These solutions share many properties with the zero-potential case, being also indexed by 
the initial value of the scalar field fy. The shape of the solutions is similar to the fi = 0 case and we again found no 
nodes in the scalar function. In Figure 1 we plotted a typical n = 4 regular solution with n = 0.3. Similar solutions 
exist also in the black hole case. In this context, we have found that the black hole solution found in jj| corresponds 
to a n = 4, k = — 1 configuration with a particular choice of (/ j , (f>h)■ 

As expected, the mass M of the solutions increases with /i while the maximal value of <f>i decreases. We will not 
address here the question of the limiting solution for /j, ^ 0, which seems to be an involved problem and a different 
metric parametrization appears to be necessary. Our preliminary results indicate that for /j / 0 the metric function 
S(ri) remains finite in this limit, while the value of M diverges. 


IV. ON THE STABILITY OF SOLUTIONS 


A. Stability of the numerical solutions 


Following the standard method, we consider spherically symmetric, linear perturbations of our equilibrium solutions, 
keeping the metric ansatz as in ©, but now the functions H, S and </> depend on t as well as r. The algebra is simplest 
if we work in the transformed frame (see section 2.3), where we have a minimally coupled scalar field, and, once the 
perturbation equations have been derived, transform back to the frame with a conformally coupled scalar field. The 
metric perturbations can be eliminated to yield a single perturbation equation for 


df = r 


("-2) 


(l-£c^ 2 ) 5 <50, 


(23) 


where 8<f> is the perturbation in the conformally coupled scalar field. For periodic perturbations ( 8cj)(t,r) = e lat 8(j){r ), 
etc), the perturbation equation takes the standard Schrodinger form 


2 T d 2 


V'F, 


(24) 


where we have defined the “tortoise” co-ordinate r* by 

dr* 

dr 


He 8 ' 


(25) 


For regular solutions, r* has values in a finite interval [0, r*i] for some r*i < oo, while for black holes, r* £ (—oo,0]. 
The perturbation potential V is given as follows (we write the formula explicitly only for the vanishing self-interaction 
potential case, for simplicity): 


V = 


He 


28 


k t ow o A 2 H 

-(n-2)(n —3)- (n _ 2)W 


Ar 2 2A£ c nr 2 
TT + (n - 2) 2 


2 n 

TT 


- (n + 2) 


/*>3 


4A ^ c n<f>(j)'r' 
Ail 


--( n -2) 3 (n-3)rV 2 


A (n - 2) 2 rV 2 

APAl 


(26) 
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where 17 is given in <UH) and A is given in 

As is often the case for AAdS solutions, care is needed in the use of boundary conditions to ensure that there is a self- 
adjoint operator in the perturbation equation El ( see, for example, the discussion in M for the Einstein-Yang-Mills 
case). We consider black hole and soliton solutions separately in this regard. 

Firstly, for black hole solutions, it is convenient to change the independent variable to y = —r* so that y £ [0, oo). 
In order to have a self-adjoint operator, we need to impose the boundary condition = 0 at y = 0, which corresponds 
to v —> oo. Secondly, for the regular soliton solutions, boundary conditions need to be imposed at both r* = 0 and 
r* = r*i as we are working on a finite interval. Suitable boundary conditions are = 0 at both the end-points, namely 
at the origin and at infinity. In both these cases, it is straightforward to check that these boundary conditions are 
sufficient to enable a self-adjoint operator to be constructed from the differential operator in m using the standard 
techniques, found, for example, in section XIII:2 of El- 

Some typical perturbation potentials El are shown in Figures 7 10. We find a complicated behaviour of the per¬ 
turbation potential depending on the values of the parameters A, <f >o and fi and the number of space-time dimensions. 
As r —■> 0, we have 


V 


(n — 4 )(n 

2^ 


1) 


0 ( 1 ), 


(27) 


so for n > 4 the potential diverges to infinity, like a standard central well potential with angular momentum (see 
Figure 8). The potential V can also be seen to vanish at the black hole event horizon (provided A > 0 and 12 > 0 
there), and, at infinity, the leading order behaviour is 


V 


r 2 n(n — 4) 

IF 


+ 0 ( 1 ), 


(28) 


so for n > 4 the potential again diverges to infinity (see Figures 8 and 10). Turning on the self-interaction potential 
V{4>) tends to increase the perturbation potential (Figure 7), as also observed in [ 5 }. 

In a limited number of cases we find that the perturbation potential El is positive everywhere (see, for example, 
some of the plots in Figure 10). In these cases, we can immediately conclude that the corresponding solutions are 
(linearly) stable, since we have a self-adjoint operator in 12411 . However, for the majority of the solutions examined, 
the potential is not positive everywhere, and in these cases we examine the zero mode solution of the perturbation 
equation (EH) , namely the time-independent solution when a 2 = 0 (see @J for further details of the zero mode method 
applied to scalar field perturbations). For all the solutions we examined, the zero inodes have no nodes (zeros). As 
the operator in El is self-adjoint, we can apply standard theorems (see, for example, section XIII:7 of [H|, or H3). 
which state that the number of nodes of the zero mode equals the number of eigenvalues a 2 of El which are less 
than zero. Therefore there are no eigenmodes which grow exponentially in time and the solutions are linearly stable. 
Our conclusion that our numerical solutions are stable is in accord with previous work ill , in which solutions 
with a scalar held having no zeros were linearly stable. 


B. Stability of the closed-form solutions 


We can also study the stability of the exact, closed form solutions on pure AdS space, given by 12211. In this case we 
keep the background AdS metric fixed and perturb simply the scalar held. The equation for the “tortoise” co-ordinate 
r* E3 can be explicitly integrated to give x = r*/(. = tan 1 ( r/t ), where x £ [0,7r/2], and the perturbation equation 
takes the form El , for periodic perturbations, with r* replaced by x, and the perturbation potential is now 

V = —- (n — 2) 2 csc 2 (a;) — - (n — 2) (n — 3) sec 2 (a;) — dp (p — 1); (29) 

where 4 p(p — 1) = (n — 2)(n + 2)/4 for the hrst type of solutions, and 4 p(p — 1) = n(n + 4)/4 for the second type of 
solutions. Since the perturbation potential is so simple in this case, we use an analytic, variational, method to study 
the stability. We dehne a functional 


= 


/ re—0 




dx 




d'b 

dx 


1 f r 5 

(d'&Y , t2 

+ / 

— TV'S 2 

0 

ll 

0 

11 

9 

\dx J 


dx, 


(30) 
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where in the second line we have integrated by parts. If we can find a test function iko such that 



(31) 


is finite, the boundary term in (|30|) vanishes and •7 r ['I'o] < 0, then there must be at least one bound state solution 
of the perturbation equation m with er 2 < 0, rendering the solutions unstable. Using \ko(a;) = sin(2a;) as our test 
function, the integral ED is finite, the required boundary term does indeed vanish, and 


•H*o] 



— 13?r 2 + 56n — 44 

u 

13?r 2 + 52n — 48 


(32) 


for the first and second type of exact solutions, respectively. Both the quadratics above are negative for all n > 4, so 
our exact, closed form solutions are unstable. This instability is perhaps a little surprising given the stability of the 
numerical solutions, however this may be understood as a result of the particular form of the (negative-) scalar field 
potential. If we perturb the scalar field slightly, there is no nearby solution of the form El for it to become. The 
nearest numerical solution will require a finite change in the scalar field at some value of and this shows up as a 
linear instability. However, it might be reasonable to expect that the scalar held would settle on a stable, non-trivial 
solution rather than radiating away to infinity. 


V. CONCLUSIONS 


In this paper we have studied the Einstein-scalar system in various space-time dimensions, with a conformally 
coupled scalar held and a negative cosmological constant. We hnd both regular and black hole solutions, generalizing 
the black hole solutions of [§). 

Both types of solutions are shown to be linearly stable, apart from some exact, discrete, closed-form solutions on 
pure AdS, which are linearly unstable. As with previous solutions Hi, this stability can be readily understood in 
terms of the Breitenlohner-Freedman bound fbll . which, in n dimensions, states that scalar helds in pure AdS are 
stable if their mass-squared satishes the inequality 0 


m BF > 


A(n 
2 {n 


JO 

2 )’ 


(33) 


noting that A < 0 so the bound is for negative mass-squared. In our case, with zero self-interaction potential, the 
“effective” mass is given by 


£ c R(r —► oc) 


An 

2 (n — 1) ’ 


(34) 


For all n > 3, it is the case that £ c I?(r —» oo) > mr BF , implying the stability of our solutions. 

There are various interesting applications of these solutions. The soliton solutions may be of interest in the 
gravitational collapse of scalar fields in AAdS. Critical collapse of a conformally coupled scalar field has been studied 
in flat space fa, but not, to date at least, in 3+1 dimensions in AAdS. Soliton solutions of the type found in this 
paper do not occur in flat space, so their presence in AAdS may change the phenomenology of gravitational collapse, 
since there are no longer just the end-point possibilities of empty space or a black hole. However, since the solitons we 
have found here are stable, they cannot be the critical solutions, unlike the situation for Einstein-Yang-Mills solitons 
in asymptotically flat space [lfj . 

It would also be of interest to calculate the mass and action of these solutions, and to look for possible applications 
within the context of the AdS/CFT correspondence. The thermodynamics of black holes with a conformally coupled 
scalar field with a positive or zero cosmological constant has already yielded some surprises (see, for example, [16j for 
a review), prompting a detailed study of the thermodynamics when the cosmological constant is negative. 

Finally, it would be interesting to extend the results derived in G2 within the isolated horizon formalism, to the 
more general case considered in this work. 

We hope to return to these questions in a future publication. 
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Figure 1 . The functions m(r), and <f>(r) are plotted as functions of radius for two typical n = 4 regular solutions. 

The scalar potential in one case is V(<f>) = n/2 </> 4 , with /r = 0.3, while V((f>) = 0 for the second solution. The central 
value of the scalar held is 4>o = 2. 
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Figure 6. The parameters M and e 2S ^ rh ' > are shown as a function of 4>{rh) for k = — 1 topological black hole 
solutions in several spacetime dimensions. 
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Figure 7. The perturbation potential m is plotted for the typical n = 4 regular solutions shown in Figure 
1 . 
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Figure 8. The perturbation potential pi is plotted for two typical n — 5, 6 regular solutions. For both the 
solutions, the scalar potential V((f>) = 0, and the central value of the scalar field is <po = 2. 
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Figure 9. The perturbation potential 123 ) is plotted for the typical n 
Figure 2. 


4 black hole solutions shown in 
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Figure 10. The perturbation potential 12fill is plotted for typical n = 5,6 black hole solutions with the value 
of the scalar field on the event horizon <po = 1.8. 













